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In this Letter, we propose a novel information-geometric characterization of chaotic (integrable) 
energy level statistics of a quantum antiferromagnetic Ising spin chain in a tilted (transverse) external 
magnetic field. Finally, we conjecture our results might find some potential physical applications in 
quantum energy level statistics. 
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I. INTRODUCTION 

o , 

Research on complexity [l[ has created a new set of ideas on how very simple systems may give rise to very complex 
behaviors. Moreover, in many cases, the "laws of complexity" have been found to hold universally, caring not at all 
for the details of the system's constituents. Chaotic behavior is a particular case of complex behavior and it will be 
: the object of the present work. 

(— I ' In this Letter we make use of the so-called Entropic Dynamics (ED) 2]. ED is a theoretical framework that arises 
O |. from the combination of inductive inference (Maximum Entropy Methods (ME), Q) and Information Geometry 
(IG) The most intriguing question being pursued in ED stems from the possibility of deriving dynamics from 
purely entropic arguments. This is clearly valuable in circumstances where microscopic dynamics may be too far 
removed from the phenomena of interest, such as in complex biological or ecological systems, or where it may just 
^ [ be unknown or perhaps even nonexistent, as in economics. It has already been shown that entropic arguments do 
account for a substantial part of the formalism of quantum mechanics, a theory that is presumably fundamental 0] ■ 
Perhaps the fundamental theories of physics are not so fundamental; they may just be consistent, objective ways 
to manipulate information. Following this line of thought, we extend the applicability of ED to temporally-complex 
' (chaotic) dynamical systems on curved statistical manifolds and identify relevant measures of chaoticity of such an 
("s^ \ information geometrodynamical approach to chaos (IGAC). 

^sO ' The layout of this Letter is as follows. In section II, we present an introduction to the main features of our IGAC. 

' In section III, using our information geometrodynamical approach to chaos and following the results provided by 
(•~^ . Prosen Q , we propose a novel information-geometric characterization of chaotic (integrable) energy level statistics of 
a quantum antiferromagnetic Ising chain in a tilted (transverse) external magnetic field. We emphasize that we have 
OO omitted technical details that will appear elsewhere. However, some applications of our IGAC to low dimensional 
Ci-i? ■ chaotic systems can be found in our previous articles 

0,i,i,[I3l- Finally, in section IV we present our final remarks. 
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II. THEORETICAL STRUCTURE OF THE IGAC 

IGAC arises as a theoretical framework to study chaos in informational geodesic flows describing, physical, bi olog ical 
or chemical systems. It is the information-geometric analogue of conventional geometrodynamical approaches 
where the classical configuration space F^; is being replaced by a statistical manifold Ms with the additional possibility 
of considering chaotic dynamics arising from non conformally flat metrics (the Jacobi metric is always conformally 
flat, instead). It is an information-geometric extension of the Jacobi geometrodynamics (the geometrization of a 
Hamiltonian system by transforming it to a geodesic flow ,13]). The reformulation of dynamics in terms of a geodesic 
problem allows the application of a wide range of well-known geometrical techniques in the investigation of the solution 
space and properties of the equation of motion. The power of the Jacobi reformulation is that all of the dynamical 
information is collected into a single geometric object in which all the available manifest symmetries are retained- the 
manifold on which geodesic flow is induced. For example, integrability of the system is connected with existence of 
Killing vectors and tensors on this manifold. The sensitive dependence of trajectories on initial conditions, which is 



'Electronic address; |carlocafaro2000@yahoo.it| 



2 



a key ingredient of chaos, can be investigated from the equation of geodesic deviation. In the Riemannian [ll| and 
Finslerian [T^ (a Finsler metric is obtained from a Riemannian metric by relaxing the requirement that the metric be 
quadratic on each tangent space) geometrodynamical approach to chaos in classical Hamiltonian systems, an active 
field of research concerns the possibility of finding a rigorous relation among the sectional curvature, the Lyapunov 
exponents, and the Kolmogorov-Sinai dynamical entropy (i. e. the sum of positive Lyapunov exponents) Using 
information-geometric methods, we have investigated in some detail the above still open research problem [8|. One 
of the goals of this Letter is that of representing an additional step forward in that research direction. 



A. General Formalism of the IGAC 



The IGAC is an application of ED to complex systems of arbitrary nature. ED is a form of information-constrained 
dynamics built on curved statistical manifolds Ms where elements of the manifold are probability distributions 
{P(X|9)} that are in a one-to-one relation with a suitable set of macroscopic statistical variables {9} that provide 
a convenient parametrization of points on Ms- The set {9} is called the parameter space Vq of the system. 

In what follows, we schematically outline the main points underlying the construction of an arbitrary form of entropic 
dynamics. First, the microstates of the system under investigation must be defined. For the sake of simplicity, we 
assume the system is characterized by an Z-dimensional microspace with microstates {X}. The main goal of an 
ED model is that of inferring "macroscopic predictions" in the absence of detailed knowledge of the microscopic 
nature of arbitrary complex systems. More explicitly, with "macroscopic prediction" we mean the knowledge of 
the probability function statistical parameters (expectation values, variances, etc.). Once the microstates have been 
defined, we have to select the relevant information about such microstates. In other words, we have to select the 
macrospace of the system. It is worthwhile mentioning that the coexistence of macroscopic dynamics with microscopic 
dynamics for a given physical (biological, chemical) system from a dynamical and statistical point of view has always 
been a very important object of investigation [isj . For example, in certain fluid systems showing Rayleigh-Benard 
convection [T6| , the macroscopic chaotic behavior (macroscopic chaos) is a manifestation of the underlying molecular 
interaction of a very large collection of molecules (microscopic or molecular chaos). Moreover, macroscopic chaos 
arising from an underlying microscopic chaotic molecular behavior has been observed in chemical reactions. In order 
to study underlying dynamics in rate equations of chemical reactions, a mesoscopic description has been adopted, 
which is given by a set of transition probabilities among chemicals. In such a description, the underlying dynamics of 
macroscopic motion is that one of stochastic processes and the evolution of probability distribution of each chemical 
is investigated jl7]. 

For the sake of the argument, we assume that the degrees of freedom {xk} of the microstates {X} are Gaussian- 
distributed. They are defined by 2Z-information constraints, for example their expectation values iif. and variances 

(xk) = fJ-k and (^(^{xk - {xk)f'j^ ^ = cr^. (1) 

In addition to information constraints, each 1-dimensional Gaussian distribution pk{xk\^J■kl ffe) of each degree of 
freedom Xk must satisfy the usual normalization conditions, 

dxkPk {xk\^J,k, (7k) = 1 (2) 

— oo 

where 

Pk {xkW, ^k) = (2^)"^ exp . (3) 

(More correctly, we should label the Gaussian probability distribution in ([3]) as P(xfe|/i^, Ufc), where P indicates 
a specific Gaussian distribution selected from a generic parametric family of distributions). Once the microstates 
have been defined and the relevant (linear or nonlinear) information constraints selected, we are left with a set of l- 

i 

dimensional vector probability distributions p (^|9) = pk {xklfJ-k^ '^fe) encoding the relevant available information 

k=l 

about the system where X is the ^-dimensional microscopic vector with components {xi,...,xi) and 9 is the 21- 
dimensional macroscopic vector with coordinates {fii,..., fii; ui,..., ai). The set {9} defines the 2/-dimensional space 
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of macrostates of the system, the statistical manifold Ms- A measure of distinguishability among macrostates is 
obtained by assigning a probability distribution P {X\Q) 3 Ms to each macrostate 9 . Assignment of a probabihty 
distribution to each state endows Ms with a metric structure. Specifically, the Fisher-Rao information metric g^i, (8) 

5^,(9) = J dXp{x\e)d^\ogpix\e)d,\ogpix\e), (4) 

with n, v — 1,..., 21 and 9^ — , defines a measure of distinguishability among macrostates on Ms- The statistical 
manifold Ms, 

Ms = |p(^|9) = J[pk{xkW, ^fe)| , (5) 

is defined as the set of probabilities {p(X|9)} described above where X £ R^^, 9 G Vq — [X^ x X„]^^ . The 
parameter space Vq (homeomorphic to Ms) is the direct product of the parameter subspaces 2^ and 2^, where 
(unless specified otherwise) = (— oo, + oo)^ and = (0, + oo)^. It is worthwhile pointing out that the possible 
chaotic behavior of the set of macrostates {9} is strictly related to the selected relevant information about the set 
of microstates {X} of the system. In other words, the assumed Gaussian characterization of the degrees of freedom 
{xk} of each microstate of the system has deep consequences on the macroscopic behavior of the system itself. 

Once Ms and Vq arc defined, the ED formalism provides the tools to explore dynamics driven on Ms hy cntropic 
arguments. Specifically, given a known initial macrostate ©(imtiai) (probability distribution), and that the system 
evolves to a final known macrostate 9'*^"^''', the possible trajectories of the system are examined in the ED approach 
using ME methods. 

We emphasize ED can be derived from a standard principle of least action (of Jacobi type) . The geodesic equations 
for the macrovariables of the Gaussian ED model are given by nonlinear second order coupled ordinary differential 
equations, 

dr2 ^ ''P dT dT ^ ' 

The geodesic equations in ^ describe a reversible dynamics whose solution is the trajectory between an initial 9^'"'*"*'' 
and a final macrostate 9^^"'*^'. The trajectory can be equally well traversed in both directions. Given the Fisher-Rao 
information metric, we can apply standard methods of Riemannian differential geometry to study the information- 
geometric structure of the manifold Ms underlying the entropic dynamics. Connection coefficients F^^, Ricci tensor 
i?^^, Riemannian curvature tensor Rf^upa, sectional curvatures ICmsj scalar curvature TZms^ Weyl anisotropy tensor 
W^vpa , Killing fields ^'^ and Jacobi fields JP can be calculated in the usual way [ll|, . 

To characterize the chaotic behavior of complex entropic dynamical systems, we are mainly concerned with the 
signs of the scalar and sectional curvatures oi Ms, the asymptotic behavior of Jacobi fields JP on Ms, the existence 
of Killing vectors ^'^ (or existence of a non-vanishing Weyl anisotropy tensor, the anisotropy of the manifold under- 
lying system dynamics plays a significant role in the mechanism of instability) and the asymptotic behavior of the 
information-geometrodynamical entropy (IGE) Sms (see It is crucial to observe that true chaos is identified by 
the occurrence of two features 'l^: 1) strong dependence on initial conditions and exponential divergence of the Jacobi 
vector field intensity, i.e., stretching of dynamical trajectories; 2) compactness of the configuration space manifold, 
i.e., folding of dynamical trajectories. The negativity of the Ricci scalar TZmsj 

TIms = R^.upag'^'g''^ = J2 (ep, e^) . (7) 

p/cr 

implies the existence of expanding directions in the configuration space manifold Ms- Indeed, since TZms the sum 
of all sectional curvatures of planes spanned by pairs of orthonormal basis elements {cp = Oq^ }, the negativity of the 
Ricci scalar is only a sufficient (not necessary) condition for local instability of geodesic flow. For this reason, the 
negativity of the scalar provides a strong criterion of local instability. Scenarios may arise where negative sectional 
curvatures are present, but the positive ones could prevail in the sum so that the Ricci scalar is non-negative despite 
the instability in the flow in those directions. Consequently, the signs of ICms ^^'e of primary significance for the 
proper characterization of chaos. 

A powerful mathematical tool to investigate the stability or instability of a geodesic flow is the Jacobi-Levi-Civita 
equation (JLC equation) for geodesic spread [ll|, 
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The JLC-equation covariantly describes how nearby geodesies locally scatter and relates the stability or instability of a 
geodesic flow with curvature properties of the ambient manifold. Finally, the asymptotic regime of diffusive evolution 
describing the possible exponential increase of average volume elements on A^^ provides another useful indicator of 
dynamical chaoticity. The exponential instability characteristic of chaos forces the system to rapidly explore large 
areas (volumes) of the statistical manifold. It is interesting to note that this asymptotic behavior appears also in the 
conventional description of quantum chaos where the entropy (von Neumann) increases linearly at a rate determined 
by the Lyapunov exponents. The linear increase of entropy as a quantum chaos criterion was introduced by Zurek 
and Paz [l8|. In my information-geometric approach a relevant quantity that can be useful to study the degree of 
instability characterizing ED models is the information geometrodynamical entropy (IGE) defined as @ , 



Sm. (t) =^^lim logV^v!^ with Vm. (^) = ^ J dr' J ^<P^Q 



(9) 



\Ms 



and g = |det (5piy)|. IGE is the asymptotic limit of the natural logarithm of the statistical weight defined on A^s and 
represents a measure of temporal complexity of chaotic dynamical systems whose dynamics is underlined by a curved 
statistical manifold. In conventional approaches to chaos, the notion of entropy is introduced, in both classical and 
quantum physics, as the missing information about the systems fine-grained state [loj . For a classical system, suppose 
that the phase space is partitioned into very fine-grained cells of uniform volume A?;, labelled by an index j. If one 
does not know which cell the system occupies, one assigns probabilities pj to the various cells; equivalently, in the 
limit of infinitesimal cells, one can use a phase-space density p{Xj) = Then, in a classical chaotic evolution, the 
asymptotic expression of the information needed to characterize a particular coarse-grained trajectory out to time r 
is given by the Shannon information entropy (measured in bits) [l9[ , 

^iiaticii^ fdXp {X) log, (p {X)Av)^~Yl Pi Pj - ICr. (10) 

where p {X) is the phase-space density and pj — ^ is the probability for the corresponding coarse-grained trajectory. 

'^dassicai^ is the missiug information about which fine-grained cell the system occupies. The quantity K. represents the 
linear rate of information increase and it is called the Kolmogorov-Sinai entropy (or metric entropy) (/C is the sum of 
positive Lyapunov exponents, /C = J2j )■ ^ quantifies the degree of classical chaos. 



III. APPLICATION OF THE IGAC TO QUANTUM ENERGY LEVEL STATISTICS 

The relevant indicators of chaoticity within the IGAC are the Ricci scalar curvature TZm,, (or, more correctly, the 
sectional curvature JCms)i the Jacobi vector field intensity Jms ^^"^ IGE once the line element on the curved 
statistical manifold Ms underlying the entropic dynamics has been specified. In what follows, selected two special 
line elements, we focus exclusively on the asymptotic temporal behavior of the IGE Sm,, related to the ED arising 
from them. 

We apply the IGAC to study the entropic dynamics on curved statistical manifolds induced by classical probability 
distributions of common use in the study of regular and chaotic quantum energy level statistics. In doing so, we suggest 
an information-geometric characterization of a special class of regular and chaotic quantum energy level statistics. 

As we said previously, we have omitted technical details that will appear elsewhere. However, our previous works 
(especially ([^j)) may be very useful references in order to clarify the following application. Recall that the theory 
of quantum chaos (quantum mechanics of systems whose classical dynamics are chaotic) is not primarily related 
to few-body physics. Indeed, in real physical systems such as many-electron atoms and heavy nuclei, the origin of 
complex behavior is the very strong interaction among many particles. To deal with such systems, a famous statistical 
approach has been developed which is based upon the Random Matrix Theory (RMT). The main idea of this approach 
is to neglect the detailed description of the motion and to treat these systems statistically bearing in mind that the 
interaction among particles is so complex and strong that generic properties are expected to emerge. Once again, this 
is exactly the philosophy underlining the ED approach to complex dynamics. It is known that the asymptotic behavior 
of computational costs and entanglement entropies of integrable and chaotic Ising spin chains are very different [Q] . 
Here Prosen considered the question of time efficiency implementing an up-to-date version of the t-DMRG for a family 
of Ising spin ^ chains in arbitrary oriented magnetic field, which undergoes a transition from integrable (transverse 
Ising) to nonintegrable chaotic regime as the magnetic field is varied. An integrable (regular) Ising chain in a general 



5 



homogeneous transverse magnetic field is defined through the Hamiltonian Tircguiar (0, 2), where 



-2 n-1 



-H {K, hy) = + {h^a^ + hVay) . (11) 

j=o j=o 

In this case, the computational cost shows a polynomial growth in time, ^j^'^^suiai) ^^-j while the entanglement 

entropy is characterized by logarithmic growth, 

^regular (0,2)= 5^on^^euma„n C log T + c' (12) 

The constant c depends exclusively on the value of the fixed transverse magnetic field intensity B± , while c' depends 
on B± and on the choice of the initial local operators of finite index used to calculate the operator space entanglement 
entropy. Instead, a quantum chaotic Ising chain in a general homogeneous tilted magnetic field is defined through 
the Hamiltonian Tichaotic (1, 1), where Ti. is defined in (jlip . In this case, the computational cost shows an exponential 
growth in time, ^j^^^^^"*"^) (f'j ^(x°° exp(/Cgr), while the entanglement entropy is characterized by linear growth, 

^chaotic (1, 1) = 5^011 Neumann ^q'''- (1^) 

The quantity JCq is a constant, asymptotically independent of the number of indexes of the initial local operators 
used to calculate the operator space entropy, that depends only on the Hamiltonian evolution and not on the details 
of the initial state observable or error measures, and can be interpreted as a kind of quantum dynamical entropy. 

It is well known the quantum description of chaos is characterized by a radical change in the statistics of quantum 
energy levels [2l|. The transition to chaos in the classical case is associated with a drastic change in the statistics of 
the nearest-neighbor spacings of quantum energy levels. In the regular regime the distribution agrees with the Poisson 
statistics while in the chaotic regime the Wigner-Dyson distribution works very well. Uncorrelated energy levels are 
characteristic of quantum systems corresponding to a classically regular motion while a level repulsion (a suppression of 
small energy level spacing) is typical for systems which are classically chaotic. A standard quantum example is provided 
by the study of energy level statistics of an Hydrogen atom in a strong magnetic field. It is known that level spacing 
distribution (LSD) is a standard indicator of quantum chaos [2^ . It displays characteristic level repulsion for strongly 
nonintegrable quantum systems, whereas for integrable systems there is no repulsion due to existence of conservation 
laws and quantum numbers. In Q, the authors calculate the LSD of the spectra of Tiicguiar (0, 2) and Tichaotic (1, !)■ 
They find that for Hrcguiar (0, 2), the nearest neighbor LSD is described by a Poisson distribution. For Hchaotic (1, 1), 
they find the nearest neighbor LSD is described by a Wigner-Dyson distribution. Therefore, they conclude that 
'^regular (0, 2) and T^chaotic (1, 1) indeed represent generic regular and quantum chaotic systems, respectively. We will 
encode the relevant information about the spin-chain in a suitable composite-probability distribution taking account 
of the quantum spin chain and the configuration of the external magnetic field in which they are immersed. 

In the ME method 3], the selection of relevant variables is made on the basis of intuition guided by experiment; it 
is essentially a matter of trial and error. The variables should include those that can be controlled or experimentally 
observed, but there are cases where others must also be considered. Our objective here is to choose the relevant 
microvariables of the system and select the relevant information concerning each one of them. In the integrable case, 
the Hamiltonian Tircguiar (0, 2) describes an antiferromagnetic Ising chain immersed in a transverse homogeneous 
magnetic field -Btransvorsc = B± with the level spacing distribution of its spectrum given by the Poisson distribution 

PPoisson (a;A|AiA) = — cxp (-— ) , (14) 



Ma V Ma/ 

where the microvariable xa is the spacing of the energy levels and the macrovariable /i^ is the average spacing. The 
chain is immersed in the transverse magnetic field which has just one component B± in the Hamiltonian Tircguiar (0, 2). 
We translate this piece of information in our IGAC formalism, coupling the probability to an exponential bath 

(exponential) / i \ • i 

Pb (a^slMs) given by 

^^exponential) ^^^j^^^ ^ ±_ / ^\ ^^^^ 

Mb V Ms/ 

where the microvariable xb is the intensity of the magnetic field and the macrovariable /i^ is the average intensity. 
More correctly, x b should be the energy arising from the interaction of the transverse magnetic field with the spin i 

particle magnetic moment, xb — —fj- ■ B — \—^B cos ip\ where ip is the tilt angle. For the sake of simplicity, let us 
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set /X = 1, then in the transverse case ip = and therefore xb = B = B±. This is our best guess and we justify it 
by noticing that the magnetic field intensity is indeed a relevant quantity in this experiment (see equation (|12p ) and 
its components (intensity) are quantities that are varied during the transitions from integrable to chaotic regimes. 
In the regular regime, we say the magnetic field intensity is set to a well-defined value (xb) — fJ-B- Furthermore, 
notice that the Exponential distribution is identified by information theory as the maximum entropy distribution if 
only one piece of information (the expectation value) is known. Finally, the chosen composite probability distribution 
p(intograbio) (^XA, xsj/iy^, /i^) cucoding relevant information about the system is given by. 



p(i„tegrable) (x^, a^sj/i^, A*b) 



1 






cxp 













(16) 



Again, we point out that our probability (jl6p is our best guess and, of course, must be consistent with numerical 
simulations and experimental data in order to have some merit. We point out that equation p6p is not fully justified 
from a theoretical point of view, a situation that occurs due to the lack of a systematic way to select the relevant 
microvariables of the system (and to choose the appropriate information about such microvariables) . Let us denote 
the two-dimensional curved statistical manifold underlying our information geometrodynamics. The line 



integrable 



M 

element ds? 



integrable 



on M 



integrable 



is given by. 



J 2 _ ,2 

""'integrable ~ "^Poisson 



ds 



1 
Mi 



Exponential — ,,2 '^I^A + ,,2 ^I^B' 



1 



(17) 



Applying our IGAC (see ([1|)) to the hue element in ([T7]), we obtain polynomial growth in ^nd logarithmic 

IGE growth, 



V 



(integrable) r^oo ^^^^^,^yc,a^ ^(integrable) r^oo ^ ^ 



M 



(18) 



The quantity c/g is a constant proportional to the number of Exponential probability distributions in the composite 
distribution used to calculate the IGE and c'jq is a constant that depends on the values assumed by the statistical 
macrovariables fi^ ^^nd fiB- Equations in (jlSp may be interpreted as the information-geometric analogue of the 
computational complexity £)("^suiar) ^^-^ entanglement entropy ^regular (0, 2) defined in standard quantum 

information theory, respectively. We cannot state they are the same since we are not fully justifying, from a theoretical 
standpoint, our choice of the composite probability ()16p . 

In the chaotic case, the Hamiltonian Tichaotic (1; 1) describes an antiferromagnetic Ising chain immersed in a tilted 
homogeneous magnetic field Stilted = B± B± + By with the level spacing distribution of its spectrum given by the 
Poisson distribution pwigner-Dyson (s^aIa^a) 



PWigner-Dyson {^a\I-^a) 



(19) 



where the microvariable x'^ is the spacing of the energy levels and the macrovariable /x^is the average spacing. The 
chain is immersed in the tilted magnetic vector field which has two components -Bx and _B|| in the Hamiltonian 
"^chaotic (1, 1)- We translate this piece of information in our IGAC formalism, coupling the probability (fTO)) to a 

^ . (Gaussian) / / i / I \ ■ i 

Gaussian (3^bIA*_Bi "b) given by. 



(Gaussian) i I \ I I \ 

Pb {.XB\^XB, (^b) = 



1 



exp 



^"b 



(20) 



where the microvariable x'g is the intensity of the magnetic field, the macrovariable fi'^ is the average intensity of the 
magnetic energy arising from the interaction of the tilted magnetic field with the spin i particle magnetic moment, 
and a'g is its covariance: during the transition from the integrable to the chaotic regime, the magnetic field intensity 
is being varied (experimentally). It is being tilted and its two components {B± and -B||) are being varied as well. 
Our best guess based on the experimental mechanism that drives the transitions between the two regimes is that 
magnetic field intensity ( actually the microvariable fiB cos ip) is Gaussian-distributed (two macrovariables) during 
this change. In the chaotic regime, we say the magnetic field intensity is set to a well-defined value (x'b) — IJ-'b with 

covariance a'g = (^{x'b — (2;^))^^ Furthermore, the Gaussian distribution is identified by information theory as 
the maximum entropy distribution if only the expectation value and the variance are known. Therefore, the chosen 
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composite probability distribution p(chaotic) {x'j^, a;^|/^^, a'g ) encoding relevant information about the system is 
given by, 



vr (2t 



-Xj^ exp 



Ms J 



(21) 



Let us denote three-dimensional curved statistical manifold underlying our information geometrody- 



(chaotic) 



namics. The line element rfSchaotic on A^^ 



"''chaotic ~ "''Wigncr-Dyson 



is given by, 

^■^Gaussian = ~72^/^A ^ 72 '^Ms ^ 72 '^'^B' 
f^A 



(22) 



Applying our IGAC (see ([!])) to the line element in (P^ . we obtain exponential growth for V^^°*"^ and linear IGE 
growth 



-y(chaotic) 



CiG exp (/C/gt) , 



(chaotic) 



(r) 



JCigt. 



(23) 



The constant C/g encodes information about the initial conditions of the statistical macrovariables parametrizing 



elements of M. 



(chaotic) 



The constant /C/g, 

3 dSM, (t) 



JC 



IG 



dT 



r— >oo 

lim 


^log( 


•^A^s (^) 


)] 


T — >00 




Jms (0) 





dot 



(24) 



is the model parameter of the chaotic system and depends on the temporal evolution of the statistical macrovariables. 
It plays the role of the standard Lyapunov exponent of a trajectory and it is, in principle, an experimentally observable 
quantity. The quantity Jms (''') Jacobi field intensity and Aj may be considered the information-geometric 

analogue of the leading Lyapunov exponent in conventional Hamiltonian systems. Given an explicit expression of 
/C/G in terms of the observables /i^ and /Xg and erg, a clear understanding of the relation between the IGE (or 



/C/g) and the entanglement entropy (or JCq) becomes the key point that deserves further study. Equations in ([23]) 

are the information-geometric analogue of the computational complexity and the entanglement entropy 

^chaotic (1, 1) defined in standard quantum information theory, respectively. This result is remarkable, but deserves a 
deeper analysis in order to be fully understood. 

One of the major limitations of our findings is the lack of a detailed account of the comparison of the theory 
with experiment. This point will be one of our primary concerns in future works. However, some considerations 
may be carried out at the present stage. The experimental observables in our theoretical models are the statistical 
macrovariables characterizing the composite probability distributions. In the integrable case, where the coupling 
between a Poisson distribution and an Exponential one is considered, and are the experimental observables. 
In the chaotic case, where the coupling between a Wigner-Dyson distribution and a Gaussian is considered, /z^ 
and //g and dg play the role of the experimental observables. We believe one way to test our theory may be 
that of determining a numerical estimate of the leading Lyapunov exponent Amax or the Lyapunov spectrum for 
the Hamiltonian systems under investigation directly from experimental data (measurement of a time series) and 
compare it to our theoretical estimate for Aj [23|. However, we are aware that it may be extremely hard to evaluate 
Lyapunov exponents numerically. Otherwise, knowing that the mean values of the positive Lyapunov exponents are 
related to the Kolmogorov-Sinai (KS) dynamical entropy, we suggest to measure the KS entropy /C directly from 
a time signal associated to a suitable combination of our experimental observables and compare it to our indirect 
theoretical estimate for /C/g from the asymptotic behaviors of our statistical macrovariables [2^. We are aware 
that the ground of our discussion is quite qualitative. However, we hope that with additional study, especially in 
clarifying the relation between the IGE and the entanglement entropy, our theoretical characterization presented in 
this Letter will find experimental support in the future. Therefore, the statement that our findings may be relevant 
to experiments verifying the existence of chaoticity and related dynamical properties on a macroscopic level in energy 
level statistics in chaotic and regular quantum spin chains is purely a conjecture at this stage. 



IV. CONCLUSIONS 



In this Letter, we proposed a theoretical information-geometric framework suitable to characterize chaotic dynamical 
behavior of complex systems on curved statistical manifolds. Specifically, an information-geometric characterization 
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of regular and chaotic quantum energy level statistics appearing in a quantum Ising spin chain in external magnetic 
field was presented. It is worthwhile emphasizing the following points: the statements that spectral correlations of 
classically integrable systems are well described by Poisson statistics and that quantum spectra of classically chaotic 
systems are universally correlated according to Wigner-Dyson statistics are conjectures, known as the BGS (Bohigas- 
Giannoni-Schmit, 25] and BTG (Berry- Tabor-Gutzwiller, j^Sj) conjectures, respectively. These two conjectures are 
very important in the study of quantum chaos, however their validity finds some exceptions. Several other cases may 
be considered. For instance, chaotic systems having a spectrum that does not obey a Wigner-Dyson distribution 
may be considered. A chaotic system can also have a spectrum following a Poisson, semi-Poisson, or other types of 
critical statistics [13]. Moreover, integrable systems having a spectrum that does not obey a Poisson distribution may 
be considered as well. For instance, the Harper model would represent such a situation. Moreover, it is worthwhile 
pointing out that not every chaotic system characterized by entropy-like quantities growing linearly in time has a 
spectrum described by a Wigner-Dyson distribution. Well-known examples presenting such a situation are the cat 
maps [28j and the famous kicked rotator [291] where its spectrum follows a Poisson distribution in cylinder represen- 
tation and a Wigner-Dyson in torus representation but the properties of entropy-like quantities are the same in both 
representations (at least classically). All these cases are not discussed in our characterization. 

Therefore, at present stage, because of the above considerations and because of the lack of experimental evidence 
in support of our theoretical construct, we can only conclude that the IGAC might find some potential applications in 
certain regular and chaotic dynamical systems and this remains only a conjecture. However, we hope that our work 
convincingly shows that this information-geometric approach may be considered a serious effort trying to provide a 
unifying criterion of chaos of both classical and quantum varieties, thus deserving further research and developments. 
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